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We re-examined Shlyakhter's analysis of the Sm data in Oklo. With a special care of minimizing 
contamination due to the inflow of the isotope after the end of the reactor activity, we confirmed that 
his result on the time-variability of the fine-structure constant, |d/a|;$10~ 17 y -1 , was basically correct. In 
addition to this upper bound, however, we obtained another result that indicates a different value of a 2 
billion years ago. We add comments on the recent result from QSO's. 
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More than 60 years ago, Dirac startled the science 
community by saying that the gravitational constant, G, 
is not a true constant, but varies as the inverse of the 
cosmic time t, and hence G/G ~ £~ 1 B ) At the present 
epoch we should expect 

| ~ -to 1 ~ -io- 10 y-\ (i) 

because the present age of the universe is to ~ 10 10 y. 
Unfortunately this number is so small that almost any 
noise would easily wipe out the real effect if there is any. 

No solid evidence for the change has ever been re- 
ported. We have heard only about upper bounds 
<10- u y -1 ^ Since these are already below the predic- 
tion ([jj), it seems as if we have lost a motivation for fur- 
ther searching for time-variability of constants. We em- 
phasize, however, that we have now another "faith," uni- 
fied theories, according to which having time-dependent 
coupling constants is a rather natural consequence, likely 
to show themselves in ways quite different from what 
Dirac expected before. 

This can be understood if we recognize that there are 
two different levels of the theories. At one of the levels we 
have a fundamental theory, like string theory as one of the 
best known candidates, while theories which apply to our 
realistic world is at another level. In the former theory 
we are supposed to have the fundamental coupling con- 
stants which are truly time- independent. On the other 
hand, the "effective" coupling constants which we mea- 
sure are derived basically from the coupling constants at 
the deeper level. But this process is accompanied with 
participation of scalar fields that are expected to vary 
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slowly as the universe evolves. This is the reason why we 
naturally expect time variation of such non-gravitational 
coupling constants like the electromagnetic charge e, or 
the fine-structure constant a = e 2 / (Airch) , as well as its 
strong interaction analogue, a s . 
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Table 1 shows a brief summary of the results of the 
past efforts for d/a and d s /a s B* a Most of them are up- 
per bounds, except for the last one from distant QSO's 
due to Webb et ai, who claim to have discovered the 
evidence of changing a for the first time.tP We will 
come back to this issue later. We notice; (i) the re- 
sults are several orders of magnitude below the level of 
to ~ 10 y , (ii) "Oklo phenomenon " yields much 
better upper bounds than others. We are going first 
to discuss the second item based on our own recent re- 
analysis. ET 

Oklo is a name of a uranium mine in Gabon, West 
Africa. This name is, however, remembered as the place 
where the evidence of natural reactofs. was discovered 
and verified for the first time in 1974.113) By natural re- 
actors we mean that self-sustained fission reactions oc- 
curred naturally some 2 billion years ago, lasting millions 
of years. Also this was a re-discovery of the theoreti- 



cal prediction due to Kuroda in 1956H Among several 
reasons for this apparently strange phenomenon to have 
occurred, most crucial is that the abundance of 235 U was 
much higher (~ 3%) in 2 billion years ago than it is 
(0.7207%). Since 1974 extensive studies have been made 
on the remnants of fission products in Oklo. 

Then in 1976, a Russian physicist A. Shlyakhter 
pointed out that the measurement of the abundance of 
the isotope 149 Sm should be useful to determine how 
much nuclear reactions in the Oklo times might have 
been different from what they are today. The natural 
abundance of 149 Sm is now 13.8%, while the measure- 
ment in the reactor zones showed much smaller values. 
For this deficit one might blame the consumption due to 
the neutron- absorption process: 

n+ 149 Sm^ 150 Sm + 7 . (2) 

This scenario is in fact supported reasonably well by the 
use of the data taken from today's laboratories. This 
implies that the process (||) 2 billion years ago should 
not have been very much different from what it is. 

The analysis can be made even more precise by noting 
that the reaction (||) is dominated by a resonance that 
lies as low as E r = 97.3meV, which is more than 7 orders 
of magnitude smaller than the typical mass scale ~ lMeV 
for most of the nuclear reactions. From this point of view, 
the occurrence of this resonance must be a consequence of 
nearly complete cancellation between the two effects; the 
repulsive Coulomb force and the attractive nuclear force. 
The former is proportional to a, while the latter could 
depend on a s probably in a more complicated manner. 
Suppose a was different from the present value by Aa. 
No matter how small Aa might be, it could still result 
in an appreciable amount of change AE r in the residual 
portion E r , particularly in terms of the fractional ratio 
AE r /E r . The same should be true for the cross section 
tTi4g of the reaction (|^) . 

This is an ideal example of an "enhancement mech- 
anism," which makes it possible to measure something 
very small. We find, as an illustration, choosing AE r = 
— lOmeV and T — 20°C gives about 10% increase of the 
cross section, where we have assumed a thermal equilib- 
rium of the neutron flux. 

Unfortunately, however, his paper was too short to tell 
any details of the data, leaving questions on how reliable 
his final result was. For this reason .we decided to redo 
what he did previously by ourselves. lP We soon realized 
that major source of errors comes from the following sit- 
uation. 

Each natural reactor came eventually to an end prob- 
ably in some million years. Even after this time, certain 
amount of isotopes might have migrated from outside 
into the core, due to the effect of weathering and other 
related phenomena. This amount, however, has nothing 
to do with the nuclear interaction we are interested in. 
We wanted to minimize this "contamination" as much as 
we could. For this purpose we used five samples taken 
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Fig. 1 Thermally averaged effective cross section for 
n+ 149 Sm ^ 15 ° Sm+7B ) The horizontal two lines rep- 
resent the range of the observed effective cross section 
(91 ±6)kb. 



from deep underground, also with the care of geologists 
expertise. In this sense we managed to have the data of 
reasonably good quality. 

We set up differential equations for the evolution of 
the system of 235 U, 147 Sm, 148 Sm, 149 Sm. By solving them 
we determined CT149 = (91 ± 6)kb for the reaction (||), 
where the "effective cross section" is defined by a = 
\J (4/7r)(T/To)<7, making it easier to handle the cross sec- 
tions obeying the 1/u-behavior with Tq = 20.4°C. 

Also we used improved methods to determine the tem- 
perature estimate; T = (200 - 400)°C. Using these to- 
gether in Fig. 1, we finally obtain two ranges of the 
allowed value of the change of the resonance energy (dif- 
ference from today's value): 

{9 ±11 meV, 
for the right-branch range, 
-97 ± 8 meV, 
for the left-branch range. 

(The result for the right-branch range replaces Eq. (27) 
in reference 9, where the erroneous lower end was shown 
inadvertently.) The presence of two solutions comes from 
the fact that the cross section is sharply peaked around 
AE r w —(50 — 100)meV, depending on the temperature, 
giving possibly two intersections with a horizontal line 
representing a value of the observed cross section. The 
right-branch range, situated to the right of the peak of 
the cross section, covers the zero, whereas the left-branch 
is away from zero. 

Damour and Dyson carefully analized the Coulomb 
energies of 149 Sm and 150 Sm, finding that Aa should 
result in AE r = (Aa/a)A4 c with the estimate A4 C w 
— l.lMeV.li3 In this way we arrive at the two ranges of 



the fractional change of a: 

Aa_ f -(0.8 ±1.0) x 10- 8 , 
a ~ \ (0.88 ±0.07) x 10~ 7 , 

and further dividing by —2 x 10 9 y also at the fractional 
rate of change: 



a 
a 



(0.4 ±0.5) x 10- 17 y-\ 
-(0.44 ±0.04) x 10~ 16 y 



(The signs in Eqs. (35) and (36) of reference 9 should 
also be reversed.) 

The result from the right-branch gives a null result 
with an upper bound, as usual. We also find that our 
result agrees well with Shlyakhter's previous conclusion, 
thus confirming his "champion" result. The agreement 
to this extent seems, however, rather accidental because 
it is unlikely that the data as good as ours was available 
in the early years of Oklo studies. 

The left-branch range indicates, as it stands, the res- 
onance energy, and hence the value of a, was indeed 
different from today's values. We tried to see if we 
could eliminate this choice by examining other isotopes, 
155 Gd, 157 Gd. The result is less conclusive due to the 
stronger effect of contamination, nevertheless indicating 
a consistency with the null result of 149 Sm. It still seems 
short of coming to a final conclusion, unless some other 
isotopes, probably 113 Cd, are included in our analysis. 

Taking contamination effect more seirously seems in- 
teresting even for 149 Sm. If we assume the contamination 
of 4%, for example, we find AE r = 2 ± 12 meV which 
replaces the previous result with contamination ignored, 
hence bringing AE r = and hence Act = more inside 
the allowed ranges corresponding to the lcr-level of the 
cross section. 

The same type jof re-analysis had been made by 
Damour and Dyson,0 who obtained the range that cov- 
ers basically both of our two ranges connected without 
separation between them. This is due to the larger errors 
in the cross section, which seems to come from contami- 
nation that affected most of the earlier data. 

We add a few comments on the recent analysis on 
QSO's.tP They improved the analysis of past studies 



based on the alkali doublets by including the effect of 
relativistic corrections to the atomic spectra, obtained 
from 49 QSO's, with the look-back time as old as 0.4-0.9 
of the age of the universe. They summarized the result 

as 

Ao 

= (-0.78 ±0.18) x 10~ 5 , 

a 

which is obviously nonzero, though some of the unknown 
systematic errors are still suspected. Comparing with 
this, our Oklo result (even the nonzero result) is about 
2 orders of magnitude smaller, but the look-back time 
is also smaller (0.16) than the aobve mentioned 0.4-0.9. 
Combining all together we would have a non-monotonic 
and even oscillating time dependence, which could be 
interesting from a theoretical point of view. 

There are suggestions that even such traditional phe- 
nomena like big-bang nucleosynthesis and anisotropy in 
cosmic microwave background allow certain room for 
time variation of alB Searching for time-variability of 
a will increasingly deserve further attention. 
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